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Abstract: There is a substantial literature on testing for the equality of the 
cumulative incidence functions associated with one specific cause in a compet- 
ing risks setting across several populations against specific or all alternatives. 
In this paper we propose an asymptotically distribution-free test when the al- 
ternative is that the incidence functions are linearly ordered, but not equal. 
The motivation stems from the fact that in many examples such a linear or- 
dering seems reasonable intuitively and is borne out generally from empirical 
observations. These tests are more powerful when the ordering is justified. 
We also provide estimators of the incidence functions under this ordering con- 
straint, derive their asymptotic properties for statistical inference purposes, 
and show improvements over the unrestricted estimators when the order re- 
striction holds. 

1. Introduction 

In the competing rislcs model, an experimental unit or subject is exposed to several 
risks, but the actual failure is attributed to exactly one of several possible distinct 
and exclusive types, called causes of failure. The available data is in the form of 
(T, 6), where T is the time of failure (from any of the causes) and S G {1,2,.. .} 
is the cause of failure. Let S denote the survival function (SF) of T. To determine 
the probability of failure by time t from a certain cause, one defines the cumulative 
incidence function (GIF), a subdistribution function (SDF), by 

(1.1) F,{t)^P[T<t,6 = i]= [ S{u)dk,{u), 

Jo 

i = 1, 2, . . . , where Aj is the cumulative cause specific hazard function corresponding 
to the ith risk. We assume that S is continuous. 

Suppose that we have k populations and in each population the units are ex- 
posed to several competing risks including one common risk whose severity we wish 
to compare across the populations. Without loss of generality, let the risk from 
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cause 1 be of interest and let the remaining causes be lumped as other causes. For 
i = 1, 2, . . . , and j = 1, 2, let Fij(t) be the cumulative incidence function due to 
cause j in population i. Let Tu be the time of failure of the Ith miit in the ith 
population and let Su be the indicator of the corresponding cause of failure. With- 
out censoring, the available data is in the form of independent random samples, 
((r,i, (Sii), (ri2, (5,2), . • . , (Ti„,, (5i„,)), i = 1,2,..., A:. If Cu is a censoring variable 
acting on the Ith unit of population i and Lu = min{Ti;, Cij}, then the data is in 
the form of {{Lu, Su)}, where Su G {0, 1, 2} with Su = corresponding to censoring 
of this observation. For censored data. Gray [7] and Pepe and Mori [14] considered 
the test of 

(1.2) H„ : Fn{t) = F2i{t) = ■ ■ ■ = Fki{t) for alH 

against all alternatives, using an integral of weighted differences of the empirical 
hazard rates and of the empirical GIF's, respectively. Note that it is not necessary 
for the other causes to be the same across the populations. Luo and TurnbuU [13] 
also developed an integral test for equality of all GIF's across the populations when 
the competing risks are the same for each population. This paper also contains 
many other related references. Lin [11] developed a 2-sided Kolmogorov-Smirnov 
(KS) type test for the case k = 2 for the same problem and provided procedures 
for simultaneous confidence intervals for Fn, i = 1,2. This test is asymptotically 
distribution-free and consistent against all alternatives. Since the asymptotic distri- 
bution of the test statistic is intractable. Lin [11] proposed a resampling technique 
for his test of (1.2) against all alternatives when k = 2. 

In many of the examples given in the literature, intuitively we would expect that 
{Fii} is linearly ordered and the empirical evidence seems to support this. In this 
paper we consider a sequential KS type test of i?o in (1-2) against Hi — Hq, where 

(1.3) i?i <F2i < ••• <Ffci. 

This test is more powerful than the test against all alternatives when the linear 
ordering holds, as is to be expected. We also provide estimators of {Fn} under the 
restriction, derive their asymptotic properties, and provide improved simultaneous 
confidence band procedures under the order restriction. 

For two linearly ordered distribution functions (DF's), Brunk et al. [2] found the 
nonparametric maximum likelihood estimator (NPMLE) whose asymptotic distri- 
bution was derived by Praestgaard and Huang [16]. This asymptotic distribution 
is too complicated for further statistical analyses. The NPMLE for the /c-sample 
case for continuous DF's is unknown. Hogg [9] had suggested a simple estimator 
that can be applied to k linearly ordered DF's. For the 2-sample case, Rojo and Ma 
[20] and Rojo [19] showed that Hogg's estimator generally has a lower mean square 
error (MSE) than the NPMLE at all quantiles. Rojo [18] and Rojo [19] showed that 
the limiting distributions of the estimators are the same as for the unrestricted case 
for k = 2 when the linear ordering is strict. El Barmi and Mukerjee [3] (hereafter 
referred to as EBM [.:>]) studied the limiting distributions for a general k, and when 
the ordering is not necessarily strict, providing procedures for asymptotic simulta- 
neous confidence bands and KS type hypothesis tests. It turns out that all of their 
results hold when the DF's are replaced by SDF's in the uncensored case. For the 
censored case, however, the asymptotic distributions of the test statistics are in- 
tractable when the DF's are replaced by SDF's. In contrast, a similar KS type test 
statistic has an asymptotic distribution that is the supremum of a rescaled Brown- 
ian Motion when testing whether all the GIF's of a single population are equal (El 
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Barmi et al. [5], El Barmi and Mukerjee [4]). We use the method in Lin [It] for 
a sequential 1-sided KS type test in the fc-sample case. The two main purposes of 
this paper are to extend Lin's test in Lin [ ' ] for the 2-sample case to the /c-sample 
case and to demonstrate how statistical inferences can be improved when the order 
restriction (1.3) holds. 

In Sections 2 and 3 we consider the uncensorcd and the censored cases, respec- 
tively. In each case we provide estimators of F21, . . . , F^i that satisfy (1.3), 
prove their consistency, derive their asymptotic properties, show the improvements 
over the empiricals in terms of asymptotic MSE (AMSE), and provide improved 
confidence bands and hypothesis testing procedures. The properties of the estima- 
tors can be proven by simply replacing the DF in EBM [-i] by an SDF. However, 
the confidence band and hypothesis testing procedures in the censored case require 
different treatments. In Section 4 we give an example of our procedures with real 
life data. In Section 5 we make some concluding remarks. 

The symbols =4>,-i, and = denote converges weakly to, converges in distribution 
to, and is equal in distribution to, respectively. 

2. The uncensored case 
2.1. The estimators 

Consider independent random samples from k life distributions, each subject being 
exposed to two causes of failure — cause 1, common to all k populations, and cause 
2, all "other causes," that need not be the same for all populations. Thus, the CIF's 
for the ith population are Fn and Fi2, given by (1.1). Let Tu be the time of failure 
of the Ith unit in the ith population and let Su e {1, 2} be the corresponding cause 
of failure. Let {Tii,Sii), {Ti2, 5i2), . . . , (Tin;, <5mi), i = 1, 2, . . . , /c, denote the samples. 
The empirical estimate of Fij is given by 

-, "i 

(2.1) F,j{t) = - V I{Ta < t, 5a = j), j = 1, 2, z = 1, 2, . . . , fc. 

Peterson [I j] shows that this is also the NPMLE that is clearly strongly uni- 
formly consistent. For a fixed t, let Fi{t) = {Fii{t), F2i{t), . . . , Fki(t))^ and let 
n — (ni, 712, ■ ■ ■ , nk)^ . The restricted estimator of Fii{t) is obtained from the iso- 
tonic regression of Fi with respect to the weight vector n: 

(2.2) i^*(0 =i?„((Ai,...,A-i)^|^),:(i), l<«<fc, 

where A = {u ^ : ui < U2 < ■ ■ ■ < Ufc}, and E^{u\A) denotes the least squares 
projection of u onto A with the weight vector w. It follows from Robertson et al. 
[17] (herafter referred to as RWD [17]) that 

(2.3) F*-^ = maxmin^Un[Fi; r, s] 

r<i s>i 

where 

s s 

Avn[Fi;r, s] = ^ niFn /Nrs, where Nrs = ^ ni. 

l—r l—r 

RWD [17] gives a comprehensive treatment of isotonic regression. It can be easily 
verified from the properties of isotonic regression that F*-^^ is nondecreasing, right 
continuous, and takes values in [0,1]. 
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2.2. Asymptotic properties 

From Corollary B, page 42 of RWD [17], we have 

max \F*^(t) - Fait)\ < max \F,iit) - Fait)\ for each t. 

l<i<k l<i<k 

The following strong uniform consistency is an easy consequence of this. 
Theorem 2.1. If Fii,F2i, ...,Fki satisfy (1.3) then 

(2.4) P[\\F*^~Fa\\^Oasn,^^, z = 1, 2, . . . , fc] = 1. 

Next, we consider the asymptotic distributions of the restricted estimators. De- 
fine 

(2.5) Z,i„. - ^A^[F^l - and Z*i„ = - i^a], i = 1, 2, . . . , fc, 

fc 

n-i, 7m = — , and assume that lim 7i„ 7^ > for all i. 

i=l 

It follows from Breslow and Crowley [1] that 

{Zllm , Z21n2 I ■ • ■ I Zklrik => (Zll, Z21, . . . , Zkl)'^ , 

a /c-variate Gaussian process with independent components, where Zn Bi{Fii) 
for all i. Here the BiS are independent standard Brownian Bridges. Let 

Ziin, = Vn[F^i - Fii] = Zii„i/V7™ and Zn = Zn/y/Ytj « = 1, 2, . . . , /c. 
Then 

(2.6) {Zllm 1 Z2ln2: ■ ■ ■ : Zkluk =^ (Zll, Z21, . . . , Zkl)^ ■ 

We first consider the convergence in distribution at a fixed point, t. Define 

(2.7) Su = {3--F,i{t)=F,i{t)}, z = l,2,...,fc. 

Note that Su is a set of consecutive integers from {1, 2, . . . , fc} with Fji{t) — Fii{t) = 
if J £ 5it, and, as n 00, 

(2.8) ^[Fji{x) - F,i{t)] ^ c», J > Su, and v^[F,i(i) - F,i(<)] j < S,t, 
where j < {>)Sit means j < (>)/ for all I ^ Su- For 1 < r < s < fc, let 

s s 

Trsn = Ijn = Nys/n, and let Trs = lim Trsn = y^77 - 

j=r j=r 

Theorem 2.2. Assume that Fii{t) < F2i{t) < ■ ■ ■ < Fki{t) for a fixed t , then 
(^nnW, ^2inW, • • • , Z*,,^{t)f 4 {Z:,{x), Z;,{x), . . . , Z,!i(x))^, 

where 

J2{r<j<s, r,seSit} 



(2.9) ^*i(^) = Vt 



max mm - 



r<i i<s P-r 
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The idea of the proof comes from the fact that the {Su : 1 < i <k} are disjoint 
"level sets" of {Fn (i)}, and the isotonic estimator of Fi in (2.3) will possibly violate 
the ordering only within each of these level sets asymptotically. The details are given 
in the proof of Theorem 2 of EBM [:!]. 

To simplify the notation, for the rest of the paper we assume that the supports 
of Fi is [0,Ti] for some < < cx) (or [0,oo)) for all i. As discussed in EBM 
after the proof of Theorem 2, weak convergence of the starred process on n jO, t^] 
requires some restrictions. We refer the reader to this paper for possible failure of 
convergence at isolated points without the following restrictions. Let 

S.,^{j:Fji forz = l,2,...,fc. 

Consider the assumption 

(2.10) inf [Fji(a;)-F,i(t)] > for aU?7 > and j > 5"^, i 1, 2, . . . , fc- 1. 

The extension of point-wise convergence in distribution to weak convergence using 
assumption (2.10) is straightforward (see proof of Theorem 4 in EBM [''>]). 

Theorem 2.3. If (1.3) and {2.10) hold. Then 

(7^ '7* 'Z* \^ — i /" 'Z* '7* '7* 

l^lln: -^2111; ■ • ■ 7 '^felnj ^ \^\\^ ^21^ ■ ■ ■ ■• ^k\) i 

where 

ry* /— ■ S{r<j<s.r,se5i} 

Z,i=vnl max nim — ^ - - . 

Note that, if Si ~ {i}, then ^ Zn under the conditions of Theorem 2.3 as 
should be the case. 



2.3. Comparison with empirical estimators 

In this section we compare some of the properties of our estimators with those of 
the unrestricted empiricals. 

From Kelly's analysis in Kelly [1(J], we have the following interesting result. 
Suppose that X = {Xi,X2, . . . , Xi)'^ has independent components, Xi ~ N{0, af), 
erf > 0, i = 1, 2. . . . , L Define the weight vector w by 




and let X* = £'w(X|^) be the isotonic regression of {Xi} with weights {1/af}, 
constrained to lie in A as in (2.2). Then 

(2.11) P{\X*\ <u) > P{\X,\ <u) forallu>0, i = 1,2,...,L 

Now suppose that Sa, as defined in (2.7), has more than one element for some t 
with < F,i{t) < 1. Then, from the form of (2.9) in Theorem 2.2, {Z*-^(t)/^ : 
j G Sit} is the least squares projection of {Zji(t) : j S Sit} with (equal) weights 
{a^^{t) : j e S^t}. Since Zj{t) ~ N{0,Fii{t)[l - F,i{t)]/j,), we have the following 
theorem. 
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Theorem 2.4. Under the conditions of Theorem 2.3, if Su contains more than 
one element for some t with < Fnit) < 1, then 

(2.12) P[\Z*,{t)\<u\>P[\Z,,{t)\<u] 

for all w > 0, j G Su. 

Two immediate consequences of this theorem are that E[Z*^(ty\^ < E[Zii{t)]'^, 
implying an improvement in AMSE, and a possible increase of confidence coeffi- 
cients for confidence bands centered around the restricted estimators as opposed 
to the empiricals with the same bandwidths. These results can be substantially 
sharpened when fc = 2. We refer the reader to Section 4.2 in EBM [o]. 

2.4- Hypothesis testing 

Here we consider testing Hq against Hi — Hq, where Hq and Hi are given by (1.2) 
and (1.3), respectively. For fc = 2 we could use the 1-sided KS test. For the fc-sample 
test, Hogg [9] suggested a sequential pairwise test as follows. We test Hoj against 
Hij — Hoj , where, for j = 2, . . . , fc — 1, 

(2.13) Hoj ■.Fn=F2i = --- = Fji and 
Hij : Fii = = • • ■ = < Fji. 

Let 

j 

i=l 

The test statistic, %i, is defined as 

Tn = max supTjn{t), where = ^/n.J Ajn[Fji - Avn[Fi; l,j - 1]]. 

'^<j<k t 

The statistic sup^ 1jn{t) is used to test Hqj against Hij — Hqi for 2 < j < k; %i is 
the maximum of these quantities, and the test rejects Hq for large values of 7n. 
Under Hq, denoting the common GIF by Fn, we have 

Tjn = ^/ Ajn[Zjia - AVj^lZia, I, j - 1]], 

where Zi„ = (Zii„, Z2i„, . . ■,ZkinY and 7„ (7i„,72„, . . ■ HknY ■ By the contin- 
uous mapping theorem, (T2„, T3„, . . . , Tfe„)^ (T2, T,, . . . , T^)"^, where 

T, = v^[Zi, - Av^[Zi; l,j - 1]], 

1 = (71,72, ■•■ ,7fe)'^, Zi = {Zii, Z21, . . . , Zki)'^ , with Zii as defined in (2.6) , 
and Aj = lim„ Aj„, that is equal to ^jTj^i/Tj with Tj = 7*- By computing 

covarianccs, it can be shown that the 7^'s are independent, and that 

Tj^B,{Fii), 2<j<k. 
Therefore Tn converges in distribution to T, where 

T ~ max sup Tj (t) . 

2<j<k t 
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and 



P{T >t) = l 



P{snpT,{t) <t,j = 2,...,k) 



k 



= 1 



\{P{&n^Bj{F^^{t))<x) 



i=2 



< 1 



(l-e-2*')fc-i ~ (fc-l)e 



the Bonferroni approximation without the independence of {Tjn} which is quite 
good unless k is fairly large. 

3. The censored case 

In this section we consider the same set-up as in Section 2 except that there may 
be censoring in addition to the competing risks. We assume that the censoring 
mechanisms are independent of the competing risks. We denote the causes of failure 
by 0, 1 and 2, where 5 = corresponds to the observation being censored. We 
assume that the censoring variable corresponding to Tij is Cij with a continuous 
SF Gi for all i and j. In this setting, we observe {(iij, <5ij)"^}, where Lij = Tij /\Cij 
and 6ij = 5ijI(Tij < Cij). Therefore {Lij, j = l,2,...,ni} is a random sample 
from Hi — 1 — SiGi, z = 1, 2, . . . , fc. Let tt^ = SiGi and let Tm = min^ n, where Ti 
is the right endpoint of the support of Fn as defined before. 

3.1. Estimators and consistency 

In order to use the martingale formulation for asymptotics, we define the uncon- 
strained estimator of Fn by 



where, for technical reasons. Si is the left continuous Kaplan-Meier estimator of Si 



where n and A are as defined in Section 2. It is well known that (||/i||o = 
suPo<t<h \hit)\) 



(3.1) 







(3.3) 




for any b < Tm- Thus, as in the uncensored case, we have 



(3.4) F(||F* - ^ as ^ CX3, i = 1, 2, . . . , k) 



1 



from the properties of isotonic regression. 
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3.2. Weak convergence 

Let Z,i„, = - Fa] and Z*„ = ^,[F,*, -Fa], 1 < ^ < fc. From Lin [11], 

(3.5) (Zii„;,Z2i„2,...,Zfci„J'^ ^ (Zii,Z2i,...,Za.i)'^ on J|[0, r^), 

i 

where {Za} are independent mean-zero Gaussian processes with the covariances 
given by (for s <t) 

Cov{Zais),Za{t)) ^ f [1 - Fa{s) - FM][1 
Jo 

+ I [Fa{s) - Fa{u)][F,^{t) 
Jo 

By independence, Cov{Zais), Zji{t)) = for i ^ j. 

3.3. Improvements by restricted estimation 

Kelly's result in Kelly [10] leading to Theorem 2.4 depends only on normality and 
independence of {Za}- Thus, Theorem 2.4 holds in the censored case also. This 
implies that the AMSE of the restricted estimators could be reduced by using the 
restricted estimators and the confidence coefficients in asymptotic confidence bands 
centered around the restricted estimators could be more than those centered around 
the unrestricted estimators using the same bandwidths. 



F^l{t) - Fa{u)] 



dAa{u) 



Fain)] 



dKi2{u) 
ttAu) 



3.4. Confidence bands and hypothesis testing 



The distributions of Z-a's are intractable. To approximate, Lin [11], based on an 
earlier paper by Lin et al. [12], develops the following martingale formulation. In 
the following, the indices z, j and I range from 1 to fc, 1 to rii, and 1 to 2, re- 
spectively. Let Y,j{t) = I{L,j > t), N,ji{t) = /(L„ < ^ I), M,ji{t) = 
%W - lo Y^J{u)dKa{u), %{t) = Y,j{t), Nuit) = Y.%i N^Jl{t), and Ma{t) = 
Miji(t). In the counting process notation, the data {(iy, %)} is represented 
by {{Yij{-), Lin [1 1] shows that, by replacing Ma by Yl'jU Viji^iji, where 
Vij/'s are independent standard normals, the conditional distribution of 



Zait) = 
(3.6) - Fait) 



[I ' Faiu)]V^jidN,,i{u) 
* V^jidN,ji{u) + Vij2dNij2{u) 



F,jV,.j2dN,^2{u) 



given {(Yiji-), Niji{-)y\ is asymptotically equivalent to the unconditional distribu- 
tion of Za{')- This clever observation allows one to simulate probabilities involving 
Za(t) by generating random samples of Viji while fixing {(Ky (•), A^^; (•))} at their 
observed values. Lin [11] considered various forms of confidence bands using the 
transformed processes 



D.a{t) = V^,g{tmFa{t)) - </>(i^,i (*))], 
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where ^ is a known function with a continuous non-zero derivative (/)' and g is 
a weight function. By the functional delta method theorem, Dnit) is asymptot- 
ically equivalent to g{t)(j)'{Fii{t))Zii{t), that can be approximated by Dii{t) = 
g{t)4>'{Fa{t))Za{t). If P[sup,^<,<,^ \D.M > fc] = «, then 

(t){F,i{t)) ± h<t<h 

is a (1 — a) confidence band for (f){Fii{-)) on [ti,t2]. This confidence band procedure 
can be carried out for each population individually. As noted above in Section 3.3, 
one could possibly improve on the coverage probabilities with the same bandwidths 
if the bands are centered at the restricted estimators. 

One could carry out the hypothesis testing in Section 2.4 with censored obser- 
vations in exactly the same way using the same notation, but using the estimators 
given by (3.1) instead of (2.1). The hmit of the process used in testing Hq^ 
against Hij — Hqj will still be given by 

(3.7) 7^-„ = yA~:[Z,i„ - Av^^ [Zi„; 1, j - 1]] ^ 7^- = [Zy - Av-y[Zi; 1, j - 1]] 

under Hqj. The distributions of the TJ 's are intractable and they are generally not 
independent in the censored case. However, we can use the resampling scheme using 
(3.6) to approximate the distribution of Zij by ^ji/^7jn to carry out the pairwise 
tests and use the Bonferroni bound given at the end of Section 2.4 to compute the 
p- value when k is not large, say less than or equal to 10. 

4. An example 

To illustrate the theoretical results, we analyze a set of mortality data provided 
by Dr. H. E. Walburg, Jr., of the Oak Ridge National Laboratory and reported by 
Hoel [8] The data were obtained from a laboratory experiment on 82 RFM strain 
male mice who were kept in a germ- free environment (group 1) and another 99 who 
were kept in a conventional laboratory environment (group 2). Each group received 
a radiation dose of 300 rads at 5-6 weeks of age. After autopsy, the causes of death 
were classified as thymic lymphoma (cause 1), reticulum cell sarcoma (cause 2), 
and other causes which we consider here to be a censoring mechanism. Hoel [8] 
writes, "Biologists believe that both of these two diseases are lethal and that they 
are independent of one another and of other causes of death." Here, we need not 
assume the independence of the causes of death, but we do assume that they are 
independent of other causes. Intuitively, we would expect that the CIFs associated 
with all risks to be smaller in the germ- free environment. Since cause 1 is considered 
to be most the lethal of all the causes, we test Hq : Fu = F21 against Hi — Hq, 
where Hi : Fu < F21, using the test statistic described in Section 3.4. The test 
statistic we use is sup(72„(t), where 7^„ is as defined in (3.7) . Since the limiting 
distribution of this test statistic is not tractable, we use the technique described 
in Section 3.4 to simulate 10,000 replicates from the distribution Tjn- The value of 
the test statistic in this case 1.11 corresponding to a p-value of 0.676. This does 
not provide any evidence for the ordering assumption, possibly because thymic 
lymphoma is extremely lethal, and the germ-free environment does not provide 
enough protection against it. However, for estimation purposes the ordering still 
seems reasonable. The empirical and the order restricted estimators are displayed 
in Figures 1 and 2, respectively. In this case, the restricted estimators are exactly 
the same. 
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Fig 1. Unrestricted estimators of the cumulative incidence functions Fn (dotted line) and F21 
(solid line). 




200 400 600 800 1000 

time 



Fig 2. Restricted estimators of the cumulative incidence functions F^^ (dotted line) and 
(solid line). 

5. Concluding remarks 

In competing risks studies, there are several tests in the hterature for equaUty of 
the cumulative incidence function due to a specific cause across several populations 
against all alternatives. Most of these are integral tests. Lin [11] proposed a 2- 
sided Kolmogorov-Smirnov type test in the case of two populations and against 
all alternatives. In some situations it is reasonable to assume that the incidence 
functions of interest in the various populations are linearly ordered. For the 2- 
sample case, the precedure in Lin [11] could be easily modified to a 1-sided test. 
This is more powerful when the ordering holds. We show how to extend this to the 
fc-sample case using a sequence of 1-sided tests. We also provide estimators of the 
incidence functions under the order restriction, derive their asymptotic properties, 
and show improvements in the inference procedures over the unrestricted estimators 
when the order restriction is obeyed, and even when it is violated by small amounts. 
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We give a real life example to illustrate our procedure. 

Acknowledgments. The authors are grateful to the Editor and an anonymous 
referee for valuable comments that have led to an improved presentation. 
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